Spin-lattice relaxation in the mixed state of YBa2Cu307_5: 
Can we see Doppler-shifted d-wave quasiparticles? 
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We present calculations of the rate of planar Cu spin lattice relaxation in the mixed state of 
YBa2Cu307_i due to (i) vortex vibrations and (ii) electron spin-flip scattering. We emphasize 
that both mechanisms give position dependent rates due to the presence of vortices, and hence the 
magnetization recovery is characterized by a distribution of rates. We conclude that relaxation by 
vortex vibrations is too slow to be a significant factor in this material. Using a semiclassical model 
of Doppler shifted d-wave quasiparticles with a linear dispersion around the nodes, our calculation 
of the relaxation rate from electron spin-flip scattering shows partial agreement with experiment. 



I. INTRODUCTION 

The spin-lattice relaxation rate is a classic probe of 
the superconducting state. The appUcation of magnetic 
fields is central to NMR, and to study spin relaxation in 
the Meissner state, great ingenuity is required to make 
measurements using a field while insuring that the re- 
laxation takes place in the absence of field. 0] However, 
in type-II materials, magnetic field enters the supercon- 
ductor naturally in the form of vortices. The cuprate 
superconductors are extreme type-II materials in which 
the internal field is relatively uniform. 

The Slichter group has made measurements of nu- 
clear spin-lattice relaxation in the vortex state of 
YBasCuaOr-a (YBCO). § A key feature of the data 
they obtain is a broad distribution of relaxation rates. A 
distribution of relaxation rates necessarily implies spatial 
nonuniformity. (Because the relaxation rates due to two 
or more mechanisms at a given site simply add, if there 
are multiple mechanisms that are all uniform, the final 
result is also uniform.) In the mixed state of a supercon- 
ductor there is intrinsic spatial nonuniformity, namely 
the vortices. 

What mechanisms of relaxation would reflect this 
nonuniformity? Possible relaxation mechanisms can be 
divided into two classes according to whether they in- 
volve an interaction with the magnetic or the electric 
quadrupolc moment of the nucleus. Given the magnetic 
nature of vortices, it seems reasonable to focus on the 
magnetic mechanisms. It is possible that vortices induce 
strains in the crystal lattice; however, any resulting vari- 
ation in relaxation by phonons would be a second order 
effect. 

Perhaps the most direct magnetic mechanism is the 
oscillation of magnetic field caused by vortex motions. 
Given the small coherence length and the consequent 
importance of thermal fiuctuations in the cuprates, this 
mechanism deserves consideration. However, the mag- 



netic mechanism which dominates relaxation in most 
conducting materials is electron spin-flip scattering: an 
electron exchanges spin with a nucleus and simultane- 
ously scatters to a new momentum state in such a way 
as to conserve energy. In a superconductor with nodes in 
its gap, this mechanism is influenced by vortices because 
the initial and final states available to the electron de- 
pend on the supercurrent velocity |^ which is a function 
of distance from the vortex cores. In addition, there may 
be electronic states associated with the cores themselves. 
In this paper we evaluate the distribution of relaxation 
arising from vortex vibrations as well as that from elec- 
tron spin-fiip scattering. 

Most previous studies of relaxation by vortex vi- 
brations have focused on the spatial average of the rate. 
We find that this mechanism is in fact strongly position 
dependent. However, in the relatively three dimensional 
compound YBCO, on which we have focused for com- 
parison with the available experiments, even the fastest 
rates we calculate from vortex vibrations are much slower 
than what is observed experimentally. 

It is the latter mechanism, electron spin-flip scattering, 
which offers the most exciting possibilities as a window 
on the electronic states in the superconducting state of 
the cuprates. It is now widely accepted that the sym- 
metry of the superconducting order parameter in YBCO 
(and many other hole doped cuprates) is d-wave. How- 
ever, there is less consensus on whether the electronic 
excitations of the superconducting state are described by 
the conventional BCS picture of weakly interacting quasi- 
particles. Microwave conductivity measurements |7j find 
that the real part of the conductivity has a Drude-like 
form with a width which drops sharply below Tc, corre- 
sponding to increasing quasiparticle lifetimes entering the 
superconducting state. However, it is not clear how to 
reconcile the temperature independent scattering rate at 
low temperatures with any known form of scattering. In 
addition, it has been argued that recent photoemission 
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measurements are inconsistent with weakly interacting 
quasiparticles. |^ 

It was first pointed out by Volovik |^ that the Doppler 
shifting of quasiparticle energies in the presence of super- 
current flow would generate a nonzero DOS at the Fermi 
surface in a superconductor with nodes in the gap. If 
the gap rises linearly near the nodes, the DOS at the 
Fermi surface will be proportional to the superfluid ve- 
locity. The supercurrent circulating around a vortex falls 
off roughly as one over the distance from the core, but is 
cut off by the spacing between vortices. Hence, Volovik 
argued that the DOS at the Fermi surface would be pro- 
portional to one over the spacing between vortices and 
therefore to the square root of the applied field. There 
will also be a contribution to the density of states from 
states associated with the vortex cores, but this is as- 
sumed to be negligible well away from the cores. 

Specific heat provides a direct measure of the DOS. 
The first experiments were done by Moler, et al, ^ and 
since then additional measurements have been made by 
other groups providing data in a wide range of 

temperatures and fields. Theory suggests that 

there should be two scaling regimes: At low tempera- 
tures and high fields, the dominant effect is the Doppler 
shifting of quasiparticles, and hence the ^/H dependence 
first predicted by Volovik is expected. However, at high 
temperatures and low fields, occupation of quasiparticle 
states by Doppler shifting is only significant near the vor- 
tex cores, and elsewhere thermally excited quasiparticles 
dominate. In this regime the field dependence of the spe- 
cific heat is expected to be linear. Much of the data fall 
in a crossover region between these two regimes. Impu- 
rities give a nonzero DOS at the Fermi surface even in 
zero field, and therefore modify these scaling arguments. 

Specific heat measurements, however, have associated 
complications, including the nontrivial process of sepa- 
rating the electronic contributions from those of the lat- 
tice and of impurities. Moreover, it has been argued 
that the field dependence of the specific heat ob- 
served in the cuprates is not due to their d-wave nature 
but rather to variation of the vortex core size indepen- 
dent of the gap symmetry. Therefore, it is beneficial to 
look for corroboration from other probes of the DOS. A 
substantial body of work on thermal conductivity is be- 
ing developed from this perspective. |l9| - [22| 

Here we study these issues from the perspective of 
NMR. Is electron spin-flip scattering the dominant (or 
at least a significant) relaxation mechanism in supercon- 
ducting YBCO? Are the electronic states well described 
in terms of Doppler-shifted d-wave quasiparticles? To 
address these questions, we have constructed a model of 
the relaxation based on Volovik's semi-classical picture, 
focusing on states near the gap nodes. At experimen- 
tally accessible fields, the vortex cores account for at 
most a few percent of the total volume of the sample. 



and we therefore do not consider here the effect of elec- 
tronic states associated specifically with the cores. To 
facilitate comparison with the available experiments, we 
have focused on the copper (spin |) nuclei in the Cu-0 
planes. 

Our results from electron spin-flip scattering are in- 
deed spatially nonuniform and similar to what is seen 
experimentally. However, the correspondence is not en- 
tirely satisfactory. The distribution of relaxation rates 
obtained is not as broad as that seen in experiments nor 
is there as much weight at the highest rates. Further- 
more, our calculations do not show as strong a field de- 
pendence as is observed. A more detailed discussion of 
these issues appears in the conclusion. 

This paper proceeds as follows: Section (|l|) defines 
the distribution of relaxation rates and the associated 



magnetization recovery. In section (HI) we outline our 
calculation of the rate of relaxation by vortex vibrations 
and describe the results of our numerical evaluation. In 



section (IV) we present our calculation of the rate of re- 
laxation by electron spin-flip scattering and discuss its 
evaluation. Finally, section (^) presents our conclusions. 



II. THE MAGNETIZATION RECOVERY FROM 
THE RATE DISTRIBUTION 

The magnetization of the whole sample as a function 
of time is given by the integral over all possible relaxation 
rates, Wi, of the product of (i) the fraction of spins which 
relax at a given rate, P(Wi), and (ii) the magnetization 
of those spins as a function of time, M{Wi,t): 



M{t)= P{Wi)M{Wi,t)dWi 



(1) 



We assume that each set of spins characterized by a given 
Wi is described by a spin temperature. Considering, for 
example, an upper satellite {—\ to — |) inversion recov- 
ery, the magnetization of each set of spins will display 
the following three exponential decay: 

M(VFi, <) = Af(W^i, oo) [1 - 2F£{Wit)] where 
1 



£{Wit) 



^-2Wit 



10 



le-emt + ±g-i2Wit , (2) 
10 10 ' ^ ' 



M(Wi, oo) is the equilibrium magnetization of the tran- 
sition, and F is the fraction of spins which arc initially in- 
verted (one ideally, but less in real experiments). Wi, or 
equivalently 1 /Ti , is defined here to be equal to one third 
the rate at which a nuclear spin flips between the states 



TO = — I and TO 



i. In principle the up and down 



transition rates are different because the populations of 
the two states are different. However, in practice, the 
thermal energy is so much larger than the nuclear Zee- 
man energy that the population difference is very small 
and need only be kept to first order. 
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In a uniform system, P{Wi) would be a i5-function. 
However, in our system the rate of relaxation, whether 
from electron spin-flip scattering or from vortex vibra- 
tions, depends on position in the vortex lattice. P{Wi) 
is obtained by adding contributions from all positions: 



P{Wi) = — dr 5{Wi- W{r)) 

J 



(3) 



B is the average field in the sample, and $o is the fiux 
quantum. Our task is therefore to calculate W{r), the 
rate of relaxation at r. 



III. RELAXATION BY VORTEX VIBRATIONS 
A. Calculating VF"" 

First we calculate the rate of relaxation due to vor- 
tex vibrations, W""" . We assume the vortices make only 
small deviations from a perfect hexagonal lattice, and we 
consider only the case in which the external field is ap- 
plied parallel to the c axis. Furthermore, we assume that 
the vortex dynamics are overdamped. Because we focus 
on YBa2Cu307_5, which is among the least strongly lay- 
ered cuprates, we use a three dimensional formalism with 
c-axis anisotropy, neglecting ab anisotropy. 

The relaxation rate due to vortex vibrations, W^"" , is 
proportional to the correlation function of the transverse 
field: § 



W^™(r) = i72x(r,c. = 7„B) 

/ + 00 
e^"*(hi(r,i).hj.(r,0)> 
-OO 



(4) 
(5) 



where is the component of the local magnetic field 
which is perpendicular to the c axis, and where () denotes 
a thermal average. 

The magnetic field is defined by 



h + V X [A • (V X h)] = $o ^ J dr, (5(3)(r - r, 



and V • h 

where A is diagonal and 

A — A — x2 



(6) 



Xab and Ac are the magnetic penetration depths in the 
ab plane and along the c axis respectively. Solving, one 
obtains for the local transverse field j2j] 

h^{p, t) = ^oJ2j dv,p{t) e-^P-"-'(*)/„^(p) (8) 



where Yiit) = {xi{z,t),yi{z,t), z) describes the position 
of vortex i at time t and 
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(9) 



We begin by calculating 

K{v,t) = (h^(r,t) 
dp 



h^(r,0)> 



dp' 



(27r)3 J (27r)3 ' 
hx(p,t)-hx(p',0)). 



=+j(P+p')-r 



(10) 



We are assuming a perfect lattice and hence the posi- 
tion dependence is completely described by the behavior 
within the first unit cell. To obtain this, we average over 
equivalent lattice sites as well as over the length of the 
vortices: 



K{p,t) 



-E- 



dzK{r, t) 



(11) 



where A^_l is the number of unit cells in a plane perpen- 
dicular to the applied field, is the length along the c 
axis, and r = r^-|-p. is the equilibrium position of vor- 
tex i and p is a two dimensional position vector within 
the first unit cell. Performing this sum and integral we 
obtain. 



K{p,t) 



A±L, ^ 



g' 



BZ 



(2^ 



(h^(g' + V^,K,t) ■ h^(g - g' - k^, -fc„ 0)) (12) 

where A± is the area of the sample perpendicular to the 
c axis, g and g' are reciprocal lattice vectors, and k is 
restricted to the first Brillouin zone. The thermal average 
in the one phonon approximation becomes 



(7) (hi(g' + k_L,fc^,i) •hi(g-g'-kj 



g-g -k 

/a/3(g' + ^±,kz)fas{g - g' - k^ , -fc^) 

+ {z X k± ■ ep){z X k_L • es){ut{k,t)ul{k,Q)) 



(13) 



where Ui(z,i) is the displacement from equilibrium of 
vortex i at height z and time t 
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r,(i) =r° + u,(z,i) 



(14) 



and ui and Ut are the longitudinal and transverse com- 
ponents of the lattice distortions 



dz Ui(z, t)e 



— zk-] 



ui(\i,t)\i± +Ui(k, t)z X ki 



(15) 



a, P,5 = x,y and repeated subscripts are summed over. 
The correlation functions of the distortions, including 
(m^) the mean square vortex displacement, are obtained 
using the vortex lattice elastic energy and the Langevin 
equation for overdamped motion of the vortices. 

If dk 

rtelasUc - 2 J^^ (27r)3 

{ei{k^,k,)\uiik)\^ 

with 

eiikj_,k:,) = cii(k_L,A:^)fc5^ 



-et(/^±,fc,)|Mt(k)|2} (16) 



C44.{k±,kz)kl 



etik±, fcz) = ceefci + C44(kj_, fc^)fc^. 



(17) 



cii, C44 and C66 are the compression, tilt and shear moduli 
of the vortex lattice respectively. The wavevector depen- 
dences of the elastic moduli are given in reference [^ . 
The equation of motion for the vortices is 



= F, 



VISCOUS 

dui 



^elastic ^ 
5'Helasti: 
5\1, 



Fr. 



(18) 



where s is the layer spacing and 77 is the vortex viscosity 
coefficient. 77 is inversely proportional to the flux flow 
resistivity. In YBCO we estimate that it is of order 8 x 
10-6 g/(cm s) for fields of order 1 Tesla at 30 Kelvin. ||2|] 
Using these. 



{ui{k,t)u*{k,Q)) = A^L, ksT 
where Ti{k±,kz) — 



j-r,(fej_,fc,)t 
^iik±,kz) 

$o£;(fc_L,fc^) 

2Bt] 



and likewise for the transverse component. Also, 
(u^) « keT 



47r$oC66C44(fcBZ, f ) 



(19) 
(20) 

(21) 



where k'^^ = gives the size of the magnetic Brillouin 
zone. Using the parameters for YBCO, which I list in 
detail below, we find (it^) is of order 10^ for fields of 
order 1 Tesla at 30 Kelvin. Combining (|l9|), ([T^), and 
( 12) , and performing the Fourier transform in time from 
( |) , we obtain 



K{p,L0) 



AB^ksT-q 



„/ JBZ 



dk 

(2^ 



g-i|g'+k^|^{u^)g-i|g'+k^-g|^(«^) 

klfapig' + k±,kz)fas{g' + k± -g,fcz) 
(kj_ • ep){k± ■ es) 



2Briuj 



h (k±, kz)f 
(z X k_L • e^)(z X k^ • e^) 



+ 



[et(k_L,fcz 



2Briuj 



(22) 



Multiplying by ^7^ and setting w = 7„i3 one has a com- 
plete expression for the position dependent rate of relax- 
ation by vortex vibrations. 



B. Discussion 

In strongly layered materials, it is possible to make re- 
liable approximations about which terms in this expres- 
sion (or its two dimensional equivalent) will dominate the 
final result. However, in the relatively three dimensional 
material YBCO, such approximations are not sufhcicntly 
accurate to yield useful results. We have therefore calcu- 
lated the required sums numerically. In doing so, we use 
lattice parameters from YBCO: a c-axis lattice constant 
of 11.68 A and (neglecting ab anisotropy) an in-plane lat- 
tice constant of 3.855 A. Furthermore, we use an ab-plane 
correlation length of 16 A, an ab-plane penetration depth 
of 1500 A, and an anisotropy ratio, F, of 5: Ac — T\ab 
and = Cab/r. Notice that the momentum sums are 
cutoff by J. Our results at 30 Kelvin for two different 
fields are shown in Figure |l|. In both cases, over most 
of the sample the rate of relaxation is extremely small. 
Near the vortices the rate peaks sharply, especially at 
lower fields. However, even at the top of these peaks the 
rate of relaxation is much slower than what is observed 
exp er iment ally. 

For vortex vibration to be an effective mechanism of 
relaxation, two conditions must be met. First, there must 
be a significant component of the field transverse to the 
direction in which it is applied (hj^). Second, the fre- 
quency spectrum of the vortex motion must show signif- 
icant weight near zero frequency (i.e. at the nuclear Zee- 
man splitting). The three dimensional vortices in YBCO 
are stiffer than the stacks of pancake vortices found in 
more strongly layered superconductors. The result is that 
in YBCO both the magnitude of the transverse field com- 
ponent and the weight at zero frequency are lower than in 
two dimensional materials. However, relaxation by this 
mechanism might still be observable were it not for the 
presence of faster mechanisms. 
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FIG. 1. The calculated rate of relaxation due to vortex 
vibrations as a function of position at (a) 8 kG and 30 Kelvin 
and (b) 80 kG and 30 K. L = (2$o/V3-B)^/^ is the spacing 
between vortices. 



IV. RELAXATION BY ELECTRON SPIN-FLIP 
SCATTERING 

A. Calculating W"^ 

We now calculate the relaxation rate due to electron 
spin-flip scattering, W'^'^^ . Again, we consider only the 
case of a field apphed perpendicular to the CuO planes. 
The nuclear transition rate is the thermal average of the 
transition rate between specific electronic, as well as nu- 
clear, states: 



(23) 



where i and / are the initial and final many-body states of 
the electron system, and indicates a thermal average. 
This rate is given by the Golden Rule: 



W, 



27r 



(m,i)(n,/) 



\{nJ\V\m,i)\''SiE^,,-E„j) (24) 



where Em^i is the total energy of the nuclear state m and 
electronic state z, and likewise for Enj. The interaction, 
V, comes from the electron-nuclear hyperfine coupling. 
Using the model of Shastry pi, Mila and Rice, 



V^^elnh^ ^Ai/+(0)5_(0) + ^Bi/+(0)5_(r„)j (25) 

where r„ are the positions of the four nearest neighbor 
coppers. The values of A± and B±^ have been estimated 



by comparison with Knight shift data. |28 2 



3.06 X 10 



'19 



erg 



(26) 



We have 



1 27r . 2 



\{f\A^S^{Q) + Y.B^S^{v,,)\i)\^5{E, - Ef)) (27) 



Note that we may neglect the nuclear Zeeman splitting in 
the energy (5-function because it is so much smaller than 
the electron's. The nuclear matrix element just gives a 
factor of three. 

In order to calculate the matrix elements of the oper- 
ator 



5_(r) = *|(r)vI/T(r), 



(28) 



we write the electron operators, in terms of Bogoli- 
ubov quasiparticle operators: 




(29) 



where is the number of nuclear sites in a plane per- 
pendicular to the c axis. The values of muq and Wka are 
obtained by solving the Bogoliubov-deGcnnes equations. 
In applying the magnetic field. Ho, we make the simpli- 
fying assumption that the spatial variation which arises 
from the formation of vortices is over sufficiently long 
length scales that at each point, r, we can simply include 
a uniform superfluid flow, Vs(r), and Zeeman splitting. 
We obtain 



Uka(l") 



,-|-i(k-|-q)T 




(30) 



where 2?iq is the center-of-mass momentum of the Cooper 
pairs (q = —kp-Vs/'^vp). The sign of the square root in 
each case is determined by the requirement that Cka be 
positive: 
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(31) 



where are the normal state excitation energies, and A 
is the superconducting gap. |Q 

We can now calculate the thermal average in equation 
(p7|). The result is a sum on wavevectors k and k' over 
terms containing one of the following three ^-functions: 
'5(ekT - ek'i): ^i^ki + fk't), or ^(ekx + ek'i)- Because the 
quasiparticle excitation energies eka must, by definition, 
be positive, the latter two ^-functions restrict the result- 
ing energy integrals to a region with a width equal to the 
nuclear Zeeman splitting. The result is vanishingly small. 
We are left, therefore, with the following expression: 

\{f\A^S^{0)+Y,B^S.{rr.)\t)\H{E,-Ef) 

n.n. 



Nl 



k k' 



[A± + 2Bj_ (cos a{kx - k'J + cos a{ky - fc^))] ^ 

/(ekT)[l-/(ek'i)]'5(ekT-ek'i) 

(ukTWk'i +«kTi'k'i)^ (32) 

To evaluate this expression, we begin by approximating 
the sums by integrals. Furthermore, because of the Fermi 
functions, the important regions in k space will be near 
the intersections of each of the four nodes with the Fermi 
surface, we convert the integral over all k space to a sum 
over the four nodes of an integral centered at an inter- 
section: 



1 \ - A± f dk 2 \ ' f 



i2n) 



(33) 



where A± is the area of the sample perpendicular to the 
c axis, and kn± and fc„|| are tangent to and normal to 
the Fermi surface at node n, respectively. Because kn± 
and fc„|| are both much less than kp, we neglect them in 
calculating the cos[a(fc — k')] factors. In this case. 



COs[Q,(^kjix kji'x )] + cos[a{kny 
= 2cos(n — n ) — 



kn'y)] 



(34) 



where n, n' — 1, 2, 3, 4 label the four nodes. 

We approximate the values of and Ak as follows: 



^k ~ hvpkw, 
Ak ~ hvik± 



and 



(35) 



where vp the Fermi velocity, and vi is the slope of the 
gap near the node. This allows us to make the following 
transformation from Cartesian to polar coordinates: 



+ A2 « ^{hvFk\\y + {hvik^y = E 



Ak 



tan-f^) = 



(36) 



The Jacobian for this transformation is E/^h^vpVi] 
Thus 



1 



k k' 



1 



(2^)V {hWvi) 

f ^d^d0 I E'dE'dO' (37) 



Finally, we write E and E' in terms of the energies of the 
initial and final states of the scattered electron as given 
by equation (^ 



e = I?„t -(- (tE 
e' = D„n +a'E' 



where cr = — 1 for < e < max[0, Unj] and a 
e > max[0, Dn^] and likewise for e'; and where 

Dniir, (j)) 

1 / TT 

= ——je^Ho + hkpVs{r) cos y4> + (2n ^1)^ 
Dn'i{r, (j)) 

1 / TT 

= +-7e^^o + hkFVs{r) COS + (2n' ~ 1)^ 



(38) 
4 for 



(39) 



In these expressions, we have made two simplifying as- 
sumptions. First, as mentioned above, fc„|| and kn± are 
much less than kp and hence |k| « |k„| = kp. Second, 
we have simplified the geometry of the system by treat- 
ing the lattice as a collection of identical cylindrical cells 
each centered on a single vortex. We neglect nonuni- 
formity along the length of the vortices and allow the 
magnitude of the supercurrent velocity to vary only as 
a function of distance from the vortex core. Due to the 
asymmetry of the superconducting gap, however, some 
angular dependence remains. 

The angles 9 and 6' now enter only in the coherence 
factor: 

(ukTMk'i +i'kT''^k'i)^ 

= i [l + aa' (sin 61 sin 61' 

+ (sgnA) (sgnA') I cos 6l| | cos e'\)] (40) 

When this is integrated over 6 and 6' one obtains i(27r)^ 
for all possible combinations of signs. 
Combining all of these points, we have 



2tt 



2x2 



EE 

n n' 

E 

<t,(t'=±1 



A, 



1 1 

c^^r) (?iV«i)'2 

+ 4:B± cos(n — n')- 

nOO 

E dE E' dE' 



{2nf 
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f{Dr,T+<jE)[l-f{Dn'l+<7'E') 

5{{Dn^+aE)-{Dn'i+<7'E')) 
e{DnT+<jE)e{Dn'i+<j'E') 



(41) 



where Q is the Heaviside function. The form of this is 
simply a sum over initial and final electron states of the 
product of (i) a factor arising from the nonlocal form of 
the interaction, (ii) the density of initial states, (iii) the 
density of final states, (iv) the probability that the initial 
state is occupied, (v) the probability that the final state 
is empty, and finally (vi) an energy conserving (5-function. 
Evaluating one of the integrals, we obtain 



1 / 7e7n?i B±a^ 



16nh V h^vpvi 

'a 



{2kBTf 



+ 4cos(n — n') — 



2kBT 2kBT 



B. Discussion 



sech^x dx (42) 



In the absence of a magnetic field, there is no Zeeman 
splitting, no Doppler shifting and all positions are iden- 
tical. The integral above is then trivial, and the rate of 
relaxation at all points is 



W""^ - 



1 / 7e7n^ B±a'^ 



h vpvi 



1.4 s" 



- + 4cos(n — ; 
at 30 Kelvin 



(43) 



where we have used a ~ 3.855 A, vp ~ 1-3 x 10'' cm/s, 
and Vi ~ vp/{14:). ||l^ The probability distribution, 
P(Wi) from (^, is just a 5- function at this value. This 
fits the NQR data very well. § 

In the presence of a magnetic field, it is no longer pos- 
sible to express the rate at an arbitrary position in closed 
form. This is due to terms of the form 



sech X dx 



(44) 



where a is not zero and b is not infinity. It is useful, 
however to look at two limiting cases. 

First, far from any vortex core, the Doppler shift is 
negligibly small relative to the thermal energy. If the 
Zeeman splitting is also negligible, the rate at these po- 
sitions is the same as in zero applied field: 



esf 



resf 



W 

mm Ha=0- 



(45) 



At 30 Kelvin, the Zeeman splitting is much less than the 
thermal energy up to fields of 10 Tesla, and may be safely 
neglected. At lower temperatures, however, the Zeeman 
splitting may become significant. As the Zeeman split- 
ting is increased from zero, the effect is initially to lower 
the relaxation rate below the zero field value. Essen- 
tially this is because, while in zero field the product of 
the initial and final densities of states has the form E'^, 
with Zeeman splitting it has the form \E — Z\ \E + Z\ = 
\E^ - Z\ where Z = {^h'y^Ho)/{2kBT). As the Zee- 
man splitting is further increased, such that the zero in 
the product of density of states moves outside the range 
of the Fermi functions, the rate rises again. The mini- 
mum rate occurs at a Zeeman splitting of roughly half 
the thermal energy and is about 20% below the zero field 
rate. The rate climbs above the zero field rate when the 
Zeeman splitting reaches roughly 85% of the thermal en- 
ergy. 

It is also possible to make a good estimate of the rate 
of relaxation near the vortices where the Doppler shift is 
much larger than the thermal energy. Here we may safely 
neglect the Zeeman splitting. Because the Doppler shift 
depends on the angle between the superflow and the node 
directions, the rate of relaxation also depends on this 
angle (although the magnitude of the variation is small, 
only about 10%). The value of this angle at which the 
rate is maximum varies between tt/A at low temperatures 
and at high temperatures, with a crossover around 20 
K. At an angle of 7r/4, the Doppler shifts associated with 
two of the nodes are +Dmax/ and with the other two 
-Dmax/V2. If Dmax/{2V2kBT) >> 1, wc may neglect 
integrals from Dmax/ (2V2fcsT') to infinity and keep only 
those from zero to Dmax/{2^/2kBT). We then obtain 



' ' m n f 



1 / leln'h B±a'^ 



16nh \ h^vpvi 
{2kBT)Dl,, (l 



(64) 



1 fA, 



8\Bi 



10 s 



,12Bi 
at 30 Kelvin 



(46) 



Here, in addition to the parameters used in the zero field 
estimate, I have assumed that the maximum Doppler 
shift is roughly given by the maximum value of the su- 
perconducting gap which is approximately 200 Kelvin. 
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FIG. 2. The calculated rate of relaxation due to electron 
spin-flip scattering as a function of position at (a) 8 kG and 
30 K and (b) at 80 kG and 30 K. L = (2$o/V3B)^/^ is the 
spacing between vortices. 

Figure |^ shows the rate of relaxation as a function of po- 
sition in one vortex unit ceU for two different values of 
apphed field. To construct these sketches and the remain- 
ing figures, we have calculated the position dependence 
of the superfluid velocity from 



,(r) cx 



E 

Q 



(z X Q) sin Q • re 

1 + A2g2 



(47) 



where the magnetic penetration depth, A, is 1600 A, the 
superconducting correlation length, ^, is 16 A, and Q is 
summed over the reciprocal lattice vectors of the vortex 
lattice. The Gaussian factor provides a smooth short dis- 
tance cutoff. We scale the magnitude such that the max- 
imum value of the Doppler shift, reached at roughly 1.5^ 
from the core, is 200 K. Inside this radius, we set the rate 
of relaxation equal to its maximum value, as if there were 
a region of normal fluid in the core. This choice may be 
controversial; however, it has little influence on the final 
form of the magnetization recovery as the cores account 
for only a small fraction of the sample area. What we 
see in Figure || is that most of the area of the sample is 
relaxed at or near the thermal relaxation rate. The rate 
then rises sharply towards the vortex cores. 
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FIG. 3. The distribution of relaxation rates due to electron 
spin-flip scattering. 

The quantitative information in Figure^ is shown more 
explicitly in Figure ^. This shows the rate distribution- 
the probability that a spin will be relaxed at a given rate- 
at both field values shown in Figure 2 as well as a lower 
temperature value for comparison. In the absence of an 
applied field, the distribution would be a (5-function at 
the thermal rate of relaxation (Eqn. ^3|). In the presence 
of an applied field, a sharp peak remains at the thermal 
rate, but the distribution is spread between this and the 
large Doppler shift limit (Eqn. ^). At low fields, near 
Hci, the high rate tale is completely negligible. At higher 
fields, more weight is shifted to higher rates; however, 
even at fields of 10 Tesla (at 30 Kelvin) one must use a 
log scale to see the shape of the distribution. Because 
the lower end of the distribution is proportional to 
while the upper end is roughly proportional to TA^^^, 
for the temperature range of interest (well below l^^ax) 
the width of the distribution increases as the temperature 
is increased. 
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. 4. The calculated magnetization recovery at 8 kG and 
and data taken at 8 kG and 83 kG, all at 30 K. 

Figure ^ shows the magnetization recovery curves corre- 
sponding to the rate distributions from Fig. ^ and com- 
pares them with experimental data from Reference 
These are upper satellite inversion recovery measure- 
ments on the planar ^'^Cu nuclei of a c-axis aligned pow- 
der of YBa2Cu307_5 with the field applied perpendicular 
to the planes. The equilibrium magnetization and the in- 
version fraction have been estimated from the raw data 
to allow presentation in this format. These electron spin- 
flip scattering results come much closer to describing the 
data than did the vortex vibration work. However, the 
data exhibit faster relaxation than our calculations pre- 
dict, a broader distribution of relaxation rates (as seen 
in their more gradual decline) and a stronger field depen- 
dence. 
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FIG. 5. The rate of relaxation, divided by temperature, as 
a function of position in the resonance line at 5, 10 and 30 K 
in 4 and 8 T fields. 

Finally, in addition to these measurements of the relax- 
ation of the full resonance line, recent experiments study 
the relaxation as a function of frequency within the line. 
1^ , ^ The resonance line in the presence of a vortex lat- 
tice has a characteristic shape, shown in Figure Ti 
data taken as a function of frequency within this line 
would allow direct observation of the position dependence 
arising from Doppler shifted quasiparticle states. Figure 
^ presents the rate of relaxation divided by temperature 
expected as a function of position within the resonance 
line at two different fields and three different tempera- 
tures. At each field there is a narrow distribution of re- 
laxation rates due to the fact that the angular variation 
of the magnetic field around a vortex is different from 
the angular variation of the relaxation rate (which is also 
temperature dependent) . 
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V. CONCLUSION 

In this paper, we have examined the two most hkely 
sources of position dependent spin-lattice relaxation in 
the mixed state of YBCO. First, assuming overdamped 
harmonic vibrations of the vortex lattice with parameters 
appropriate to YBCO, we find that relaxation by vortex 
vibrations, although strongly position dependent, is too 
slow to play a significant role in the observed magnetiza- 
tion recovery. Second, we have studied electron spin- flip 
scattering within a framework of noninteracting quasi- 
particles with energies Doppler shifted by the supercur- 
rents circulating around vortices. Again our results are 
strongly position dependent. In addition, the rates we 
calculate are similar to those seen in experiments. How- 
ever, significant discrepancies remain between our calcu- 
lations and what is observed: our calculations predict 
slower rates, a narrower distribution, and less field de- 
pendence than are seen in the experiments. 

Given the emphasis we have placed on the nonunifor- 
mity of relaxation rates, it might appear that nuclear 
spin diffusion should also be considered. This is espe- 
cially so given the strong indirect nuclear dipole coupling 
through electron spins in the cuprates. If spin diffusion 
were to take place, it would smooth out the spatial varia- 
tion in relaxation rates and hence widen the gap between 
our results and the experiments. However, spin diffusion 
is very strongly suppressed in the nonuniform magnetic 
field of the vortex lattice, especially so near the vortex 
cores where the field gradient (and also the Ti variation) 
is greatest. 

The fact that this theory which involves a number of 
simplifying assumptions and experimental data with at- 
tendant complications do not agree completely is not so 
surprising. Regarding the theory, the fastest relaxation 
rates come from regions near the vortex cores. It is pre- 
cisely in this region where our quasiclassical approxima- 
tion and our assumption of linear gap variation break 
down. Hence, our theory is least reliable in its descrip- 
tion of the fastest rates. Recent numerical work may 
help clarify the behavior in this region. Regarding the ex- 
periments, impurities, inhomogeneities and practical dif- 
ficulties may influence the data: Fluctuating moments of 
magnetic impurities might increase the rate of relaxation 
of nearby nuclear sites. An inhomogeneous superfluid 
density due to defects in the grains could complicate the 
position dependence of the superfluid velocity and possi- 
bly increase the area of large Doppler shifting. Failure to 
invert spins across the full vortex-lattice field spectrum, 
particularly tricky in the presence of a distribution of 
quadrupolar splittings and grain misalignments, might 
also influence the data; although the effect in this case 
would be to narrow the observed rate distribution. 

Doing experiments on clean single crystals would of 
course settle some of these issues. However, returning to 



the issue of our approximations, what experiments would 
focus most clearly on the features which the present cal- 
culation captures? The goal is to work at a temperature 
and field which maximize the percentage of spins, V, for 
which the magnitude of the Doppler shift is (i) larger than 
the thermal energy (and the Zeeman splitting) but (ii) 
much less than the maximum gap. The first condition is 
necessary for the effects of the Doppler shifting to be ob- 
servable. The second condition ensures that our approx- 
imations are valid. Clearly going to low temperatures 
helps to satisfy the first condition. This is limited by the 
fact that the Zeeman splitting will eventually become 
larger than the thermal energy. However, for tempera- 
tures above 1 Kelvin, the fields which maximize V are too 
small for the Zeeman splitting to be an issue. Practical 
issues, such as minimizing impurity effects, will probably 
place greater restrictions on the choice of temperature. 
The following results are based on an approximate form 
of the superfluid velocity in which Vs(r) oc (l/r) — {1/ L) 
where r is the distance from the vortex core. If we set the 
upper bound on the Doppler shift at Amax/5, at 2 Kelvin 
a field of 2 Tesla will maximize V at 70% (with roughly 
10% of the spins outside the range of our approximations 
and 20% relaxing at the thermal rate). At 5 Kelvin, a 
field of 4 Tesla maximizes V at 55% (with roughly 20% 
of the spins outside the range of our approximations and 
25% relaxing at the thermal rate). 

Finding clear confirmation of the presence of Doppler 
shifted quasiparticles would be a significant achievement 
in itself. In addition, Ti measurements might provide 
insight into the nature of the quasiparticle interactions. 
In a parallel calculation of the electron spin-flip scatter- 
ing nuclear spin-lattice relaxation rate within the frame- 
work of antiferromagnetic spin fluctuations interacting 
via short-lived quasiparticles, |3^ ] many of the qualita- 
tive features of the relaxation are similar to the nonin- 
teracting quasiparticle model presented here. However, 
some features unique to each model exist, and if observed 
could clarify which model is most appropriate to YBCO. 
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